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1 Introduction 



Experimental information on the shape of the strangeness changing scalar /o(t) and vector 
form factors in the low-energy region can be obtained from the study of fC^s-decays. /o(t) and 
f+{t) indeed enter the differential decay rates of these semi-leptonic processes. In the expres- 
sion of these decay rates /o(t) is multiplied by a kinematic factor (m^/mA-)^ with and niK 
the lepton and the kaon mass, respectively. This factor being of the order 10^^ for the elec- 
tron, only the muon mode is, in fact, sensitive to the scalar form factor which is thus harder to 
determine. Different collaborations, namely ISTRA [IJ, KLOE EUS], KTeV mH and NA48 
IIHQ have measured these kaon decays. In the analysis of their data, they parametrize the two 
form factors in terms of some free parameters. The actual number of parameters which can 
be determined from a fit to the data are, due to strong correlations between them, at most two 
for the vector form factor and one for the scalar one. Thus quadratic-, pole- and more recently 
parametrizations based on conformal mapping (denoted z-parametrization) have been used for 
f+{t) while, up to recently /o(t) was described in terms of a linear and a pole one |. While 
the pole parametrization gives comparable results for in the different experiments, the 

situation for /o was much more confused. The slope of the normalized scalar form factor, /o(t), 
varied typically between 9 ■ 10^'^ and 15 ■ lO^'^ depending on the experiments [8J. Note that 
the slope determined in this way can only be an upper limit of the true mathematical slope 
dfo{t) / dt\t=o as calculated for example within chiral perturbation theory O [T01|. Therefore 
it seemed appropriate, in particular in the case of the scalar form factor, to develop another 
parametrization which is as model independent as possible, involves only one parameter and 
determines the higher order terms in the series' expansion on physical ground. An accurate dis- 
persive representation has been constructed in Ref. [ ITI which fulfills all these properties. Two 
classes of parametrizations can hence be distinguished depending on whether or not physical 
information is used. In the first class one finds for example the pole one for the vector form 
factor and this dispersive parametrization for the scalar form factor and in the second one the 
linear and the quadratic parametrizations. Strictly speaking the z-parametrization enters also 
this latter class. However, in Ref. [fT2ll it was shown that under certain conditions it is possible 
to impose a bound on the sum of the expansion coefficients based on unitarity and the total rate 
of r — i> KtxVt- 

While the main aim of these Ki^^ experiments was to extract the CKM matrix element Vus 
the dispersive parametrization of /o provides another test of the Standard Model (SM) through 
the measurement of the only unknown parameter In C, with C the value of the scalar form 
factor at the Callan Treiman (CT) point. This was in fact the original idea behind writing such a 
dispersive parametrization. It has first been used by the NA48 collaboration leading to a rather 
small slope for /o and a 4.5cr deviation to the SM [7 |. Unfortunately the situation is still unclear 
for the scalar form factor, recent determinations of In C from KLOE [O and KTeV |fT3l lead 
to no/slight discrepancy with the SM. Due to the importance of such a measurement it is very 
important to discuss in more detail the robustness of this dispersive parametrization. This is the 
main aim of this paper together with the investigation of a similar parametrization for the vector 
form factor improving on the pole parametrization. 

^Note that in the case of the scalar form factor the lightest resonance, the k*, is rather broad. In addition, there 
is a second resonance not far away. Thus contrary to the vector form factor case, the parametrization of the scalar 
form factor in terms of one real pole has no physical motivation. 
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After introducing basic notations and properties in section |2] we will briefly review in sec- 
tion |3] the dispersive parametrization for the scalar form factor. We will then describe an analo- 
gous parametrization for the vector form factor. In section |4] we will investigate the robustness 
of these parametrizations for both, the scalar and the vector form factors. We will discuss the 
uncertainties due to the input parameters as well as the expected size of isospin breaking effects. 
In the dispersive parametrization the standard, since most likely, hypothesis of the absence of 
zeros in the form factor has been made. In section 14.31 we will question this hypothesis and 
in particular, we will study the presence of possible real or complex zeros in the form factors 
and the impact of these zeros on the value of In C. Even if the likelihood of such a scenario is 
small its study is required by the particular importance of an unambiguous test of the SM. The 
possibility of discarding zeros in the form factors from some properties of the scalar form factor 
as well as from a comparison with high energy data from r-decays will be discussed in detail. 
We will conclude in section [5] and finally present some useful expressions to simplify the use of 
the dispersive parametrization in the data analysis in the appendix. 



2 Basic definitions and properties 



The hadronic matrix element describing f^^s-decays is written in terms of two form factors 

/^-(t)and/^-(t), 

{7,{p^)\s^^u\K{pk)) = (p.+Pk), + (PK-Pn), /"""W , (2-1) 

where t = (pk — P-k)'^- The vector form factor f^'^{t) represents the P-wave projection of the 
crossed channel matrix element {Q\s'y^u\KTx), whereas the S-wave projection is described by 
the scalar form factor defined as 

m = + -T-^f-'^ii)- (2-2) 

In the following discussion we will consider the normalized form factors 

m = ^ and hit) = MI with /o(0) = /+(0) = 1 , (2.3) 

and try to describe their shape as precisely as possible in the physical region of Ki^-decays, 
mj < t < to = [niK — JtItt)^, with the pion mass. It is shown in Fig. [H together with the 
right hand cut from Ktx scattering which starts at tx-K = {^k + '^tt)^ as well as the CT point 
Aktt = — TTT'l whose value is about twice as large as tg- This point is of special interest in 
the case of the scalar form factor. Indeed, the Callan-Treiman low-energy theorem [14] predicts 
its value in the SU (2) chiral limit (where the quark masses m„ ^ vanish) at that particular point. 
As we will see, this is of great importance in testing the SM. For physical quark masses, one 
has 

C . /„(A,,) = ^J^^ + Ac. , (2.4) 

where Fk+,tt+ are the charged kaon and pion decay constants, respectively, and Act is a cor- 
rection of O {mu,d/4:7rFT,) U- It has been estimated within Chiral Perturbation Theory (ChPT) at 



^Note that for practical purposes we have introduced some isospin breaking effects already in the first term on 
the right-hand side (RHS) of Eq. ( IZ4l i. 
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Figure 1: Different energy scales involved in the analysis of the scalar and vector form factors: 
the physical region of Ke^-decays lies between the lepton mass squared mj and to = {rriK — 
ttIt^Y, ^k-k = ^1: — nT''i denotes the CT point and the right-hand cut from Kit scattering starts 
at txn = {jriK + ^tt)^- For the numerical values, mK+ and m^o have been used. 



next to leading order (NLO) in the isospin limit with the result 

A^|;0 = (_3.5±8) X 10^^ (2.5) 

The error is a conservative estimate assuming typical size corrections of O {rriu^) and O {mg) 
[fT6H for the higher orders. It should certainly hold for the neutral kaon decays which we are 
mainly interested in at present. Indeed, no large corrections to this estimate are expected due 
to the absence of n^-r] mixing in the final state which could lead to small energy denomina- 
tors. Eq. (|2.5I) is consistent with the values obtained recently, see Refs [fTOl . [fTTl and [[TSl . The 
quantities entering the expression of C, Eq. (|2.4I) . are in principle completely determined within 
QCD. Except for Act, they are studied within lattice QCD, see Ref. [ 19| for a recent overview 
of the situation. However, the actual most precise determinations of these quantities are ob- 
tained from semi-leptonic decays. Consequently, they depend on the assumption made for the 
electroweak couplings of quarks. Assuming the SM couplings, one can extract the quantity 



Bexp 



IKdl , (2.6) 



u, i/r^^(o)Ksi 

using experimental information on the ratio Tr^+ /T+ , the decay — > n^eue EOl and 

0+ 0^ transitions in nuclei fT2^. These determine respectively the first ratio, the second one 
and \ Vud\ in Eq. (|2.6I) with high precision and lead to 

= 1.2418 ±0.0039 . (2.7) 

Interestingly the knowledge of Vus is unnecessary for determining B^xp- Eq. (12.41) then becomes 

In C I SM = In Bexp + Act/ Bexp 

= 0.2166 ± 0.0034 + (-0.0035 ± 0.0080)/(1.2418 ± 0.0039) 

= 0.2138 ±0.0073, (2.8) 

where in the last line all the errors have been added in quadrature. Note that the error on In CI^m 
is rather small. Hence a precise measurement of In C in neutral i^^a-decays should allow to test 
the SM electroweak couplings by comparing the obtained value with the one from Eq. (12.81) . 
However, the Callan-Treiman point is unreachable by a direct measurement of this decay, its 
value being much larger than the end point value of the physical region, see figure [Tl For this 
reason a dispersive representation of the scalar form factor written in terms of In C as the only 
free parameter has been introduced in Ref. [llj. 
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3 Dispersive parametrization 



Let us first briefly review this dispersive parametrization of the scalar form factor based on an 
Omnes representation ll23l . see also Ref. [|24ll for an early application. We will then introduce 
a similar representation to accurately describe the vector form factor in the low energy region. 



3.1 Scalar form factor 

The dispersive representation of the scalar form factor introduced in Ref. [fTTIl follows previous 
attempts to determine the form factor in the physical region. In Ref. [|25l a coupled channel ap- 
proach was used with only one subtraction in the dispersion relation. The main aim of Ref. f TTIl 
was to stay as model independent as possible. Therefore a second subtraction has been made in 
order to minimize the bad knowledge of the high energy region in the dispersive integral. Using 
the two points t = (where by definition, /o(t) = 1) and the Callan-Treiman point A^^ as 
subtraction points leads to: 



/o(t) = exp (InC-G(t)) 



(3.1) 



with 



^ Ax.(Ak-. - t) r ds Ms) ^ (3,2) 

JtK^ s {s- AKn){s-t-ie) 



Here, (f)o{s) is the phase of /o(s). In writing Eq. (13.11) . it has been assumed that /o(t) has no 
zeros. We will come back to this point in Sec. 14. 3[ In what follows, G{t) is decomposed as 

G{t) = GkMs, t) + G,,(As, t) ± 5G{t) , (3.3) 

where the first term corresponds to an integration from the threshold tx-K up to a cut-off Ag 
which characterizes the end of the elastic region while in the second term the integration runs 
from As to oo. The choice of the value of A5 will be discussed later (see section |4?T1) . In the 
elastic or low-energy region {tx-K < s < A5) the phase is identified with the s-wave, 1=1/2 
Ktx scattering phase, 5q, according to Watson's theorem [26]. In the analysis of Ref. [fTT| , 
5q (with its uncertainty) has been taken from Ref. [27 J. There a matching of the solution of 
the Roy-Steiner equations with the K-k Kn, nn KK and tttt ^ tttt scattering data 
available at higher energies has been performed. The phase obtained in this way is in very 
good agreement with the work of Ref. [25|, see below. In the inelastic or high-energy region 
(s > A5), the phase is almost unknown. Perturbative QCD indicates [|29l that fo{t) vanishes as 
0{\/t) for large negative t. Thus from Eq. (|3.2I) . we conclude that the phase must approach 
TT asymptotically. In Ref. lITTI . the phase has been taken constant and equal to its asymptotic 
value of TT for s > A5 and an uncertainty of ±7r has been assumed. Note that this uncertainty 
is a rather conservative estimate leading to a large band going from to 27r. However, due to 
the two subtractions G{t) converges rapidly, hence G{t) is almost insensitive to the high energy 
behaviour of the phase, cf. Sec. 14.11 and the large uncertainty on the phase at high energy turns 
into a small uncertainty on G{t). For example, G'(O), which has the largest error, is given for 
A5 = 2.77 GeV2 by 

G(0) = 0.0398 ± 0.0018 ± 0.0036 ± 0.0017 , (3.4) 
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where the first/second error correspond to the error on GxJGas, respectively, and the third error 
comes form the study of isospin breaking. We will come back to the uncertainties on G{t) in 
detail in section |431 

The only free parameter, In C, in Eq. (13.11) could, in principle, be determined from the sum 
rule 



dictated by the asymptotic behaviour of /o(t), cf. Ref. f29\. However, this sum rule, which 
exhibits one less subtraction than G{t), Eq. (|3.2I) . is not precise enough to allow to determine 
In G with a good accuracy without adding any information on the high energy behaviour of the 
phase of the form factor. We will come back to the discussion of this sum rule later. Thus, 
In C is a free parameter which can be determined from experiment by fitting the i^^s -decay 
distribution with the dispersive formula for /o(t), Eq. (13.11) . For more details on the dispersive 
representation of the scalar form factor, see Ref. [[TTl . 



3.2 Vector form factor 

In connection with the recent precise measurements of the differential spectrum of the r 
Ut-Ktc decay by Belle [30] and BaBar [31 J, theoretical work has been devoted in the last few 
years to the description of /+ (t) . Form factors have been obtained in the framework of reso- 
nance chiral theory with additional constraints from dispersion relations in Refs. [|32l [33l . In 
Ref. [l34ll . a coupled channel analysis has been performed taking into account, through analytic- 
ity requirements, the experimental information on elastic and inelastic Ktt scattering from the 
LASS collaboration. All these studies impose constraints from short distance QCD as well as 
the value of the vector form factor at zero momentum transfer. A fit to the r data allows them 
to determine completely the shape of the form factor and thus to deduce a value for the slope 
and the curvature which compares reasonably well with the recent Ki^, experiments. In these 
works, the emphasis is put on the energy region of the r-decay and they are thus best suited 
for it. Here our aim is somewhat different: we want to have a very precise parametrization of 
the form factor at low energy improving on the pole parametrization usually assumed in the 
ii'fs-analysis 



which expresses the vector form factor completely in terms of a resonance described as a dis- 
crete pole dX \/t = My. This parametrization is physically motivated by the dominance of the 
K*{SQ2) resonance in the vector channel. We will add to our knowledge of the presence of 
this resonance, the properties of analyticity and a proper behaviour of the phase at threshold. 
Contrary to the analysis discussed previously, we will not be able to determine the slope of the 
vector form factor but it will be a free parameter to be determined from a fit to the /T^s-data. 
This will allow us, using a twice subtracted dispersion relation as in the scalar case, to minimize 
the effect of the high energy region in the dispersive integral over the phase of the form factor. 

In the case of the vector form factor, the value of at t = is known, see Eq. (12.31) . 
but there is no equivalence of the low-energy theorem of Callan and Treiman, Eq. (|2.4I) . Thus 




(3.5) 



hit) 




(3.6) 
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a dispersion relation for ln/+(t), this time twice-subtracted at zero, will be written. Defining 
\df+(t) /dt\t=o = A+/m^ and assuming again that the form factor has no zero, this will be 
discussed in section |431 one has: 

hit) = exp (A+ + H{t))] , where H{t) = <L T ^-Alf) , (3.7) 



ml 



TT /^^, [s — t — ie) 



with 01 (s) being the phase of At sufficiently low energies, in the elastic region, is 

given by the p-wave, 1=1/2 Kir scattering phase, ^xtJ ('^) = ^li^)- ^ detailed partial- wave 
analysis of Kn Ktx scattering in the energy range sq = (0.825 GeV)^ < s < (2.5 GeV)^ 
has been performed in Ref. ^TT\ based on high statistics production experiments. In order to 
reliably evaluate the dispersion integral, Eq. (13.71) . an accurate extrapolation of the scattering 
phase down to threshold is needed. Contrary to the s-wave case, the Roy-Steiner equations are 
not really useful for providing such an extrapolation due to the lack of relevant experimental 
results. However, the well known method due to Gounaris and Sakurai [|36ll can be used to 
directly construct a partial wave amplitude which is unitary, has the correct threshold behaviour, 
the correct analyticity properties (neglecting the left-hand cut) and reproduces the position and 
width of the K*{SQ2) as given by the PDG [21 J. It is most suited to use the inverse amplitude 
method. Defining the function D{s) via the T-matrix as 



allows to determine its discontinuity 



3 

ImD(s) = -2^^ . (3.9) 
Hence, one can write a dispersive representation for 1/T leading to 

D{s) = -2qj^^{s)- j= + P[s) (3.10) 

with P{s) being a subtraction polynomial. In the previous equations, we have used the standard 
notations, 

( . ^2 ((s - {ttik + m^Y) (g - {niK - m^ f)f^^ 

S={PK+Pn), qKn{s) = , (3.11) 

where qx-K is the absolute value of the three-momentum in the Kir center of mass frame. Thus 
the Kn scattering phase is given by: 

ctg(5i(.)) = qUs)Hs) + ^ , (3.12) 

with 

his) = -^P r . (3.13) 



./tj^^ X x/X X- S 



A minimal choice will be done here for P(s) namely 

P{s)=a + bs, (3.14) 
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since it already gives a very good description of the phase in the vicinity of the resonance 
K* (892) and down to threshold. We will consider the impact of a higher order polynomial in 
section lT.2.21 The constants a and b are determined from the mass and the width of the K* (892) 
characterized as 



-Ml, 



and 



ds 



=M2 

if* 



(3.15) 



Note that there exists in the literature another definition of mass and width in terms of the 
position of the pole in the complex plane. The latter is process independent. The uncertainties 
coming from the inputs used for Mk* and Tk* will be discussed in section 14.2.21 Another 
possibility would be to determine a and b from a direct fit to the data [35] . This would, however, 
lead to a function H(t) lying within the error bars discussed below. We checked that the phase 
constructed in this way with no free parameters leads to values of the p-wave scattering length, 
oi = 0.0183m^, agreeing with other determinations ETl [37l[38l . 



3tt 



27r- 



our model 
a..=-7 10"-' 



TT 




Figure 2: Comparison of our model for the phase of the vector form factor, Eq. d3.i2l) . with the 
coupled channel analysis ofRef ^34^ . The grey band corresponds to the assumption that above 



Ay the phase equals Tcll'^, see text. 



In p-wave scattering, inelasticity effects which imply (t)i{s) ^ 5i(s) become important at 
lower energies than in the scalar case, the mass of the vector resonance fC*(1414) being an 
indication of the start of the inelasticity. At high energy, following the same arguments on 
the asymptotic behaviour as for the scalar case, the phase will reach its asymptotic value, vr. 
Therefore, similarly to what has been done for G'(t), the function Hit), Eq. (13.71) . is decomposed 
into two parts: 

Hit) = HKA^v,t) + HasiAv,t) ± SHit) , 



with 



and 



ml t 



Slit') 



t'^it' -t) 



dt' + 



ml t 



TT 



so 



Se.p{t') 
t'^it' -t) 



dt' 



ml t 



TT 



m: 



dt' = ^ In 1 - 



t 



mt 

TV 



(3.16) 



(3.17) 



(3.18) 



J Ay t'\t' -t) t "~V Xv- 

In these equations. Ay denotes the end of the elastic region. In what follows, we will use 
Ay = (1.414)^ GeV^ and we will discuss other values for Ay in section 14.2.21 In Eq. (13.171) . 
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the analytic formula for the phase 5i, Eq. (13.121) . is used below sq. Above sq and below Ay 
the experimental points of Aston et al. B5II are used to determine 5exp and to evaluate the 
corresponding contribution to H{t) together with its uncertainties. As in the scalar case, the 
asymptotic contribution, Eq. (13.181 ). gives only a very small contribution to H(t) due to the two 
subtractions. For the error, we will take a somewhat larger band than in the scalar case. Indeed 
it was found in Ref. [|34l that a fit to the r data reproducing the rate quoted by the PDG 
[|2TI and compatible with asymptotic QCD, leads to a phase of the vector form factor reaching 
Stt at infinity. Thus we will take 



corresponding to the assumption that above Ay, 0i(s) = irl^"^. This recipe almost certainly 
overestimates the real uncertainty. The other sources of uncertainty will be discussed in section 
14.2.21 Let us give here only as an indication the value of H dXt^, point in the physical region of 
iT^s-decay where it has the largest error. For Ky = (1.414 GeV)^, one has 



where the first uncertainty is the one on Hk-k and the second one the one on Has, Eq. (|3.16l) . 
Note that as for G{t), even if the inelasticity sets in at lower energies, the uncertainty on the 
value of Has{Av,t) related to our poor knowledge of 0i(s) for s > Ay is small due to the 
two subtractions performed this time at zero. The resulting phase, 01 (t), with its uncertainties 
(grey band) is displayed in Fig. fusing the charged K* mass, Mk* = 891.66 MeV and width 
Tk* = 50.8 MeV from PDG [|2TI . It is compared with the phase obtained in Ref. [34] from two 
fits characterised by a different value of the parameter om, see Ref. [34] for more details. 

Finally, let us note that a sum rule comparable to Eq. (13.51) is implied by the asymptotic 
behaviour of the form factor. 



4 Discussion of errors and assumptions 

The aim of the dispersive representation is to describe the shape of the form factor in the phy- 
sical region of i^'^a-decays as precisely as possible and consequently to test the Standard Model 
through the value of the scalar form factor at the Callan-Treiman point. In order to achieve this, 
it is mandatory to have under control all the assumptions entering the construction of the form 
factors as well as the determination of the errors. In the following, we will discuss the choice 
of the cut-off, the dependence on the input parameters and the absence of zeros. 



4.1 Cut-off dependence 

As just discussed two energy regions are distinguished in the dispersive analysis of the Kn form 
factors, the elastic one at low energy which is very well under control and the inelastic one at 



5Has{Av,t) 



+2HasiAv,t) 
-Has{Av,t) 



(3.19) 



H{to) = (2.16 ± 0.04+°:|) X 10"^ 



(3.20) 




(3.21) 
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higher energies which is of less importance in the description of the form factors due to the two 
subtractions in the dispersive representation. As explained before, in the latter region a rough 
estimate of the phase is made, using (f)o{t) = ir ±n for the scalar form factor and = n^l'^ 
for the vector one. Let us first discuss this approximation in a bit more detail. 

In Fig. [3l we show the phase of the scalar form factor as obtained from a once subtracted 
dispersion relation by Jamin et al., Ref. [25J, as well as the phase of the Kn amplitude obtained 
by Biittiker et al., Ref. [|27l . As can be seen, in the elastic region both phases agree as they 

400 



300 

CD 

I 200 

to 



100 



"0 2 4 6 8 

s [GeV]^ 

Figure 3: Comparison of the Kit scattering phases from the amplitude and the S -matrix in the 
scalar channel extracted in Ref. 1 .27] and the phase of the form factor (preferred fit 6.10K2) 
obtained via a multi-channel analysis / l25l/ . 

should. In the region where the inelasticity sets in both phases decrease, even rather abruptly 
in the case of the phase of the amplitude, and then start to grow again. This behaviour is well 
understood and has been explained in Refs. [|39l ffOj in the case of the scalar tttt form factor. 
Even though the central value of the phase discussed in section [STTI does not have this property, 
the large uncertainty assigned to it takes it into account. Unless the form factor has a zero at 
some higher energy, as it will be discussed in the next section, no other sharp drop of the phase 
is expected and the phase will just in some way go to its asymptotic value vr at very large t as 
typically does the phase obtained by Jamin et al. in Ref. [i25il . Thus (pQ^t) = n ± n certainly 
encompasses the physical phase of the scalar form factor. 

Another source of uncertainty comes from the fact that the energies and Av where the 
inelasticity cannot be neglected any more are not very well known. For the s-wave. As was 
chosen in Ref. [1 1 J as the energy where the phase of the amplitude is experimentally found to be 
different from the phase of the ^-matrix, namely at s = (1.66 GeV)^ as shown in Fig.[3l In the 
case of the vector form factor, the _ft'*(1414) resonance can be seen as an indication of the end 
of the elastic region. We want to investigate here how sensitive G{t) and H(t) are to variations 
of the cutoffs A5 and Ay within reasonable bounds. We will concentrate on G{t), an analogous 
study on H{t) leads to similar conclusions. In Fig. HI the bands represent the possible values of 
G(0), Eq. (13.21) . and of 00), Eq. (13.51 ). as a function of the cut-off where G{0) has the largest 
uncertainty in the whole physical region of -fC^s-decays. For the s-wave a reasonable range of 



f.f. 125J 

ttK amp. 1271 

ttK S-matrix 127] 
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Figure 4: Variations ofG{0) (LHS) and G{—oo) (RHS) as a function of the end point of the 
elastic region As. Note that the scale differs in the two figures. 



values for As is (1.43 GeV)^ < A5 < (1.66 GeV)^ where the lowest value is determined 
by the iir*(1430) resonance. Within these limits Gmm(O) varies between 0.0331 and 0.0354 
and Gmax{0) between 0.0474 and 0.0442 while G'™„(-oo) varies between 0.1227 and 0.1495 
and Gmax{—oo) between 0.3599 and 0.3234 where Gmax/min{t) = G{t) ± 5G{t) and the plus 
sign corresponds to the maximum value. Hence although each part of the integral naturally 
depends on the exact cut-off value, the sum is clearly less sensitive on the exact position of 
As. We obtain, for instance, G{0) = 0.0336 + 0.0067 = 0.0403 for A5 = (1.43 GeV)^ and 
G{0) = 0.0362+0.0036 = 0.0398 for A5 = (1.66 GeV)^. As expected G(- 00), which involves 
one less subtraction, depends much more on the exact value of the cut-off and has a much larger 
uncertainty. However a comparison of the theoretical result for In G, Eq. (|3.5I) . together with its 
experimental determinations (71131 [131, which lie between 0.14 to 0.21 show that the uncertainty 
on G{t) is certainly overestimated. As already stated, a precise determination of In G from the 
sum rule, Eq. (13.51) . is in the actual state of the art not possible. 



4.2 Variation of the input parameters 

Within this section we will discuss the influence of the choice of the input parameters on the 
phase at low energies and therefore on the dispersive representation of the form factors. 



4.2.1 Scalar form factor 



At low energy the / = 1/2, s-wave Kn scattering phase, 60 from Ref. [27J is used in the 
evaluation of the function G{t). As already discussed in Ref. [fTT|, there are two sources of 
uncertainties on Sq. The first one comes from the propagation of the errors of the experimental 
inputs into the solution of the Roy-Steiner equations. The other one is due to the choice of the 
point where one matches the Roy-Steiner equations with the data. This leads to the conservative 
estimate (cf. Ref. LHJ), 

SGkA^s, t) < 0.05 GkA^s, t) . (4.1) 
Furthermore, the analysis of Ref. ||27l is performed under the assumption of perfect isospin 
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Figure 5: The function G{t) in the physical region. Each curve corresponds to one of the 
uncertainty discussed in the text added in quadrature to the previous one, in the order shown in 
the legend. The last curve gives the total uncertainty on G(t). 

symmetry. In Ref. (4V\, it has been shown that for vrvr scattering the shift in the phase due to 
isospin breaking corrections can amount up to 0.015. For the nir system, these leading effects 
are of purely electromagnetic origin. In Refs. [421 l43l it has been argued that due to a partial 
cancellation between the strong and the electromagnetic effects, the total isospin breaking effect 
for Kn should be smaller than for tttt. Awaiting a more detailed quantitative analysis [|44l|. we 
will assume here as a conservative estimate of isospin breaking effects a constant shift of the 
phase by 0.02 which is even larger than the maximal value obtained for tttt [41 J. The corre- 
sponding error on GKTri^s, is, in fact, comparable to SGk-k, Eq. (14.11) . Replacing the average 
values of m^- and throughout the solution of Roy-Steiner equations with their physical ones 
has only a negligible effect on GK^i-^s, 

The function G{t) with the different uncertainties discussed here is plotted in Fig. [51 In 
section [ATI the simplified expression for G{t) discussed in Ref. LIU is given together with the 
different errors on the parameters. 

4.2.2 Vector form factor 

To evaluate Hktt we considered two domains as explained in section 13.21 From threshold tx-K 
to So the main input parameters are the mass and the width of the K*(892). Varying them 
within the error bars given by the PDG [ }2T| has only a negligible influence on Hktt, the errors 
are at most a few 10""^ times Hktt- Another source of uncertainty in this energy region lies in 
the choice of the polynomial P(s), Eq. (13.101) . In order to estimate the effect of higher order 
terms, we added a contribution c ■ to the r.h.s of Eq. (13.141) and varied the parameter c within 
-0.15 GeV"^ < c < 0.62 GeV"^ Within this range we still obtain values for the scattering 
length, ai = 0.013 — 0.020, in agreement with other determinations Il27l l37l l38l and the 
phases agree reasonably well with the Aston data as displayed in Fig. |6l The induced error on 
Hk-k is in this case a few 10^^ times its value. Between sq and the cutoff Ay the uncertainty 
on the phase is given by the error bars of the Aston data. To determine the error induced 
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Figure 6: Kn scattering phase in the vector channel for different values of the parameter c. 
Here sq is fixed to the value (0.825 GeV)^. For comparison the data / l35l/ are shown, too. 

by the choice of the cutoff values, we impose a rather large variation for sq and Ay, namely 
(0.825 GeV)2 < sq < (1.1 GeV)^ and (1.2 GeV)^ < Ay < (1-6 GeV)^. 

The function H{t) with the different uncertainties discussed here is plotted in Fig.|7l Except 
for Has, we have symmetrized the errors which were not symmetric around the central value, 
using as error the largest value to be on the conservative side. In section \A2] a simplified ex- 
pression for H{t) is given together with the different sources of uncertainties on the parameters. 



4.3 Discussion of possible zeros 

Writing the dispersive relations, Eqs. (|3.1[|3.7|) . we have assumed that the form factors have no 
zero. This is in fact quite standard when studying form factors. Indeed, in the space-like region, 
a form factor represents the Fourier transform of a charge density. It is argued for instance 
in Ref. Il45l for the case of the electromagnetic form factor of the pion that the properties of 
the pion charge distribution should be similar to the one of the electron in the ground state of 
the hydrogen atom. This charge density being proportional to the square of the wave function, 
the corresponding form factor is positive in the space-like region. The possibility that the form 
factors have one or several zeros can, however, not be completely overruled and has been studied 
in the literature for example for the pion form factor in Ref. B6l . Let us therefore discuss it 
here for the scalar form factor. In the vector case the experimental situation is much better as 
we have seen in the introduction. Slope and curvature can be measured rather precisely. Zeros 
would essentially modify the relation between the slope and the curvature in the physical region 
of /^^a-decays compared to what we have been discussing here. Considering the precision of 
the experiment, allowed zeros cannot alter much the results in the physical region. We will 
however briefly mention what happens for the vector form factor in Sec. 14.3.41 in relation with 
the T KtiUt decay. 
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Figure 7: The function H{t). Each curve corresponds to one of the uncertainty discussed in the 
text added in quadrature to the previous one, in the order shown in the legend. The last curve 
gives the total uncertainty on H{t). 



4.3.1 Real zeros 



We will consider first one real zero, Tq, to simplify. In that case the function /o(t) /(I — t/To) 
has no zero and one can proceed as in the previous section. One can thus write 

m = (i - 4) 

In order for /o(t) to behave like 1/t for large negative t l|29l . the phase of the form factor has 
now to go to 2n at infinity since one has 

lim fo{t) = const. . (4.3) 

t— > — oo 

Now there are two possibilities. The zero can be in the time-like or in the space-like region. 



InC - In 1 



G{t) 



(4.2) 



• Time-like region In the time-like region a zero on the real axis would correspond to a jump 
of the phase by vr. In fact, Eq. (14.21) can be rewritten as 

^ Juk 5 {s - AKn){s - t - ie) J \ 

Following the discussion of section l4~n this zero should occur at a four momentum larger than 
~ (1.7 GeV)^ since below, the phase of /o(t) is known and there is no indication of a zero for 
the form factor. At these energies there is a right-hand cut on the real axis and this scenario 
looks thus rather improbable. Furthermore since the region above As gives a small contribution 
to G{t) in the .fr^s-physical region, one does not expect much change in the form factor from 
the presence of such a zero. This is illustrated in Fig. [8] where a comparison is made between 
the form factor obtained from Eq. (13.11) and the one from Eq. (14.21) for Tq = 2.5 GeV^. This 
is typically the smallest possible value for a real zero in this time-like region and thus the one 
which is expected to affect most our representation. 



/o(t) = exp 



A 



K-n 
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Figure 8: Scalar form factor: without zeros plain black curve, with a zero at Tq = 2.5 GeV^, 
red dashed curve and with two complex conjugate zeros at Zq = (2.5 ± i 10) GeV^, blue 
dashed-dotted curve and at Zq = (0.1 ± i 2) GeV^ pink dotted curve. 



• Space-like region In the space-like region, there is no left-hand cut and thus a real zero 
looks more plausible even though it goes against the argument of positivity of the form factor in 
the space-like region given previously. A completely different behaviour of the form factor as 
compared with the case without zero could result if the zero was close enough to the physical 
region. This can be seen in Fig. |9] where the scalar form factor is shown for two different values 
of To, namely —0.1 GeV^ and —1 GeV^. In the first case, the slope is rather large and there 
is a maximum at a rather small t value. As Tq decreases, the slope becomes smaller and the 
maximum moves towards larger t and eventually approaches the curve without zero. This is in 
fact what happens for Tq = — 1 GeV^. The K^^^a-data, however, seem to exclude the large slopes 
we observe when the zero becomes close to the physical region, see Fig. |9l In that case, the 
high energy behaviour of the form factor is completely different, too. This has an impact for 
instance in r-decays as discussed below. 



4.3.2 Complex zeros 



Zeros could occur in the complex plane, which would in fact mean, due to the property of real- 
analyticity of the form-factor, the presence of a zero, Zq, and its complex conjugate Zq. The 
form-factor would take the following form 

(4.5) 

In order for /o(t) to behave as 1/t for large negative t {29], the phase should now go asymp- 
totically to Stt, see Eq. (14.31) . As can be seen in Fig. [8l the presence of complex zeros could 
lead to lower values for the form factor at the end of the physical region. In the case of the blue 
dot-dashed curve which corresponds to zeros with rather large imaginary parts, namely ±10 
GeV^ and a rather large real part, too, thus not close to the physical region of _ft"^3-decay, the 
difference to the curve without zero is very small as expected. The effect becomes, however. 
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Figure 9: Scalar form factor in presence of zeros located in the space-like region. Black plain 
curve: no zero, blue dashed/red dotted-dashed curves: real negative zero at Tq = — 0.1 GeV^ 
and —1 GeV^, respectively. 

much stronger when one takes a small real part as well as a smaller imaginary part, as shown 
by the pink dotted curve in Fig. [8] which corresponds to zeros at (0.1 ± i 2) GeV^. The slope 
becomes smaller and the curvature larger. Note that the same value for In C has been used in 
the determination of the curves, namely InC = 0.2138, Eq. (12.81) (the trends discussed here are 
independent of the precise value of In C) such that the three curves meet at the CT point. 



4.3.3 Properties of the scalar form factor and their constraints on the existence of possi- 
ble zeros 

It is clear from the previous discussion that the presence of zeros in the scalar form factor with 
values close to the physical region would affect the determination of In C from experiment. 
Unfortunately, it is not possible to test such a scenario with the f^^s-data since due to strong 
correlations and the sensitivity of the data it is only possible to determine one parameter for 
the scalar form factor from the fits. However, the scalar form factor possesses some properties 
which might allow to set up constraints on the presence or the absence of zeros. The different 
behaviour for t away from the physical region of i^^^a-decays could have consequences, too, for 
example in r-decays. 



• Sum rules The first properties which can be used are sum rules which can be derived from 
the knowledge of the behaviour of the form factor for t — oo. In the absence of zeros, one 
obtains the sum-rule, Eq. (13.51) . It is modified in the presence of zeros and for instance, with 
one real zero, it reads 

G(-cx)) = In C - In (^1 - . (4.6) 

The second logarithm on the right-hand side is negative so that the sum rule now leads to a 
value G{—oo) larger than InC. For example, with Tq = 2.5 GeV^ as used above, one obtains 
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ln(l — Aktt/Tq) = —0.096. Thus a good knowledge of these two quantities could provide an 
information on the presence of zeros in the form factor. In the case of complex zeros, one can 
derive two additional sum-rules namely 

° ds, /"/f^ , , = , (4.7) 

^ {l-s/Zo){l-s/Zo) 



tK 

OO 



ds- =0. (4.8) 

(l-s/Zo)(l-s/Zo) 

All these sum rules have to be satisfied simultaneously. Unfortunately, they are rather sensitive 
to the phase at high energies which, as has already been stressed, is badly known. Consequently, 
they do not provide any constraint on the presence or absence of zeros as long as one does not 
have a better knowledge of the behaviour of this high energy phase. At present, it is possible to 
find a model for the phase which satisfies simultaneously the three sum rules. 



• Soft kaon analog of the Callan Treiman theorem Another important property which pro- 
vides us with a more severe constraint is the soft kaon analog of the Callan Treiman theo- 
rem [|47l . One has 

f,(ml-mj,) = ^ + AcT. (4.9) 

A one loop calculation of the SU (3) correction in the isospin limit BUl gives = 0.03 
which is larger than its soft pion analog Act, see Eq. (12.51) . by a factor m\/m\. It is rather 
small for a first order SU (3) ® SU (3) breaking effect, which is expected to be of the order of 
25%. Note one interesting point: at NLO within the minimal not-quite decoupling electroweak 
low-energy effective theory (LEET) [ |48| , there appear in the light quark sector essentially two 
combinations of parameters of spurionic origin describing the couplings of quarks to the W- 
boson to be determined from experiment P9l [50l . While the knowledge of the scalar form 
factor at the CT point measures one combination, its knowledge at — m\ measures the 
other one. A precise determination of Kqt would thus help to settle the issue of the presence 
of right-handed couplings of quarks to the VT-boson. At present one can only give an estimate 
of the higher order contribution to Act- The expected size is of the order of ten percent. To get 
an idea of the actual size, let us look at the two-loop calculation [ISTl . where one finds: 

F Ifi 

Act = 2 - ^ - -f - -^(2Ci2 + C,,)ml{ml - ml) + A{ml - m],) + A(0) . (4.10) 

The quantities A{t) and A(0) are discussed in Ref. f5\\, and C12 and C34 are two low energy 
constants (LECs). The same combination of these LECs in Eq. (14.101) appears in the two loop 
calculation of the slope, Aq, of the scalar form factor. Taking Aq between 0.009 and 0.016 which 
encompasses the values obtained by the NA48, KTeV and KLOE collaboration and using the 
value 1.22 for Fk/Ft, as in Ref. [3T| (for an actual status on the value of F^/Ft,, see for 
example Ref. [8]) one obtains —4 x 10~^ < 2 C12 + C34 < 8 x 10"^ in agreement with the 
estimates found in the literature [[521 [531 [H [HI. This leads to -0.035 < Act < 0.11 within 
the expected order of magnitude. Assuming the Standard Model electroweak couplings, one 

^Note that in Ref. ifSTl the parametrization of A(i) is in principle only valid in the physical region. We will 
however use it here in order to get an estimate for Act- 



16 



has Ft,/ (Fr-/+(0)) = 0.8752 ± 0.0020, such that a conservative result for the normalized scalar 
form factor fo{ml — m\) is 

0.8</o(m2-m^)<l. (4.11) 

This bound constrains the allowed region of possible zeros in the scalar form factor. Real zeros 
in the space-like region which do have an impact on the form factor in the physical region are 
excluded since they drop very fast with t for t negative and violate the bound. For complex 
zeros the situation is less obvious. It is easy to calculate which are the imaginary parts the 
complex zeros should have in order that the scalar form factor at t = — m|- lies within 
the bound, Eq. (|4.1 II) . One finds that the positive imaginary part increases when the bound 
decreases. Typically for a complex zero with a real part between -0.2 and 0.2 the imaginary 
part varies between 0.7 for fo{m1 — m\) = 1.03 and 6 for fo{m1 — m\) = 0.87. Thus, with 
our present knowledge of Act, one cannot totally eliminate the presence of complex zeros very 
close to the physical region which would affect the determination of In C. 



4.3.4 High energy behaviour and the decay r Kt^Vt 

The energy distribution in the decay r KtxUt has been measured first by the ALEPH colla- 
boration [|54l , then by OPAL [|55l , and recently the Belle collaboration [30] and the BaBar one 
[|3TI presented results for different isospin combinations. The cross section is described by 

^^Yi£t^,,,^M,.l^\ (4.12, 



d^t 1287r3 -^""^ ' V K 

with the kinematic variable qxtt defined in Eq. (13.1 II) . It thus involves the scalar and the vector 
form factors. Unfortunately due to the presence of the i^'*(892) resonance which dominates 
the cross section between ~ 0.8 and ~ 1.2 GeV, one can only hope to obtain information on 
the scalar form factor from this decay channel very close to threshold. As discussed above, the 
only zeros which could lead to a clear difference in the energy distribution very close to the 
threshold are the real space-like zeros close to the i^'^js-physical region. As is shown in Fig.iTOl 
the curve is shifted towards larger values if there is no zero in the scalar form factor. The Belle 
data, seem to favor the latter even though the error bars are rather large for the lowest poinlI§. 
This is consistent with the soft kaon analog of the CT theorem which, as we have seen, excludes 
the presence of such zeros in the scalar form factor. 

On the right-hand side of the figure, the differential decay width, Eq. (14.121) . is shown in a 
broader energy range and assuming the presence or absence of zeros in the vector form factor. 
The resonance region, where the contribution of the vector form factor dominates, becomes 
more sensitive to the phase in the inelastic region such that relatively large uncertainties for the 
decay width are expected in this resonance region. The dispersive representation without zero, 
for A+ = 0.02450, and 0i = vr beyond the cut off Ay, as discussed in section \32\ describes 
amazingly well the precise Belle data in this region, see the central plain black curve on the 



We do not show the BaBar data since they are not pubUcly available. The same trends hold, however, for 
them, too. 
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Figure 10: Differential decay width r KnUr. On the left-hand side the calculations with one 
real zero in the scalar form factor close to the physical region (Tq = —0.1 GeV^j and without 
zero are compared. The vector form factor has no zero. On the right-hand side the calculation 
is done without zero and with complex zeros at Zq = (0.1 ± 12) GeV^ in the vector form 
factor Different cases have been considered corresponding to different ansdtzefor the phase as 
explained in the text. The scalar form factor has no zero. In these figures InC = 0.2138, Eq. 
dZSj) and A+ = 0.02450 as obtained from the pole parametrization with the K*(892) mass. For 
comparison the Belle data for K^t[~i't are displayed, too. 

RHS of Fig. [lOl Much larger or smaller phases at the beginning of the inelastic region seem 
to be excluded as shown on the same figure where instead of tt, for the lower plain black 
curve and Svr for the upper plain black curve have been used for 01. In the presence of complex 
zeros, on the contrary, the phase should be very small at the beginning of the inelastic region 
in order to reproduce the data, as illustrated in Fig. [lOl There the result with a phase equal to 
3 TT in the energy range from 1 .4 to 3 GeV (pink dotted curve) is compared with the one with 
a vanishing phase in that same energy range (blue dashed curve). The opposite is true for real 
space-like zeros. As the value of Tq increases, the resonance peak gets more and more washed 
out and eventually disappears for a given phase 0i . Rather highly improbable large values of the 
phase in the inelastic region becomes necessary to counterbalance the effect of the zeros. Thus 
from our study, due to the lack of knowledge of 0i in the inelastic (high-energy) region, we 
cannot completely rule out the presence of zeros in the vector form factor even though such a 
scenario does not seem very probable. We can, however, conclude as expected that zeros in the 
vector form factor which would not be totally excluded from the analysis of tau decays would 
not affect the low energy region of the vector form factor and consequently the results of the 
analysis of i^'^a-decays. 



5 Final remarks and Conclusion 



In this paper we have discussed the robustness of a precise and convenient dispersive represen- 
tation of the scalar and vector Ktx form factors. In Fig{TT]we show the scalar (left panel), Eq. 
(13.11) , and the vector form factors (right panel), Eq. (13.71 ). with all the uncertainties discussed 
above under the usual assumption of no zeros in the form factors and for two different values of 



18 




Figure 11: Dispersive representation of the scalar (left panel) and vector (right panel) form 
factor over a wide range of energies. For comparison two different values o/lnC (right) and 
A+ (left) have been used: \nC\sM is from Eq. dZgP . \nC\^AAfi = 0.1438, A^^^^ = 0.0233 are 
the central values of the NA48 experimental results [7] and A^°'^ = 0.02450 is from the pole 
parametrization with the K*{892) mass. The band takes care of all the uncertainties discussed 
in the text. 

their respective free parameters InC and A+. The figure nicely illustrates the fact that our dis- 
persive representation describes to a very high accuracy the form factor shapes in the physical 
region of ii'^s-decays. Eqs. (|3.1[|3.7| ) thus represent a very useful tool for an optimal analysis of 
the i^^s-data. As already pointed out, it allows to determine the shape of these form factors in 
an unambiguous way, contrary to other parametrizations used in the data analysis. Furthermore, 
as emphasized in Ref. [ITTI . a measurement of In C, with C the value of the normalized scalar 
form factor at the Callan-Treiman point, allows to test the Standard Model. A departure of the 
measured value from Eq. (|2.8I) would signal, under the hypothesis of no zeros in the form factor, 
a failure of the SM, as for example the presence of a direct coupling of right-handed quarks to 
W [fTTTl . We have, however, to moderate slightly the conclusion drawn there. We have indeed 
seen that the shape of the scalar form factor could be slightly modified in the highly improbable 
case where it would have zeros in a very small domain of the complex plane within or close to 
the physical region. Even though the likelihood of this scenario is very small we have not 
been able at present to totally eliminate it. Note, however, that for the vector form factor, zeros 
that would affect the dispersive parametrization in the low energy region we are interested in 
here, are excluded by the tau data. 
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A Some useful expressions 



A.l Scalar form factor 

In order to facilitate and accelerate the numerical evaluation of the scalar form factor for exam- 
ple in the experimental analysis of i^^g-decays, it is convenient to have a parametrization of 
the function G{t). As already discussed in Ref. IfTTll . within the physical region it can be very 
accurately parametrized as 



Gp{t) =xD + {1- x)d + x(l - x)k 



(A.l) 



where x = t/to, d = G{0), D = G{to). The value of k can be obtained from the constraint 
G{AKn) = 0. In Table[T]we give the central values for the parameter d, D, k together with the 
corresponding errors arising from the different sources of errors as discussed in the previous 
section. 





Central value 


5Ga 


5 Gas 


SGktt 


^G isospin 


total error 


d 


0.0398 


0.0005 


0.0036 


0.0018 


0.0017 


0.0044 


D 


0.0209 


0.0002 


0.0016 


0.0010 


0.0010 


0.0021 


k 


0.0045 


0.0000 


0.0001 


0.0002 


0.0003 


0.0004 



Table 1: Coefficients arising in the parametrization Gp with their uncertainties. 



For practical purposes and for comparison with the traditionally often used linear and quadratic 
approximations to the form factor, it is useful to list the first coefficients of the Taylor expansion 
of the form factor. They read 



J_dfo 
ml dt 

1 rf7o 



mi dt'^ 



1 d^fo 
mt dt^ 



t=o 



t=o 



mz 



^(InC-G(O)) 



mt 



^(InC- 0.0398(44)) 



m: 



m 



m. 



dfo 
dt 

2 ^ 
dt 

2 dfo 



dt 



t=o 



t=0 



t=0 



t=0 



t=0 



(4.16 ±0.56) X 10" 



771 

3^G"(0) 

^Ktt 



3 (4.16 ±0.56) X 10-4 ( ml ^ 

dt 



t=o 



+ (2.72 ±0.21) X 10"^ . 

(A.2) 



The Taylor expansion up to third order of the scalar form factor allows to well reproduce the 
exact dispersive representation in the physical region. The maximal error is 3%. 



20 



A.2 Vector form factor 



It is equally very convenient to have a parametrisation of the function H{t) in order to avoid the 
evaluation of the dispersive integral for the vector form factor. Within the physical region, the 
function H{t) can be very accurately parametrized as 



Hp{t) = Hix + H2X' 



(A.3) 



The numerical values of the parameters Hi and H2 are given in Table [2] together with the corres- 
ponding uncertainties. As discussed in the text, the uncertainties coming from the uncertainties 
on the mass and the width of the K*, 5i/Af . and 6Hr^, are completely negligible. Here again. 





Central value 


Aston 


5Hs, 


5H, 


5Ha^ 


SH,s 


total error 


Hi X 10^ 


1.92 


0.04 


0.01 


0.01 


0.06 


-0.31 


+0.63 
-0.32 


H2 X 10^ 


2.63 


0.04 


0.01 


0.02 


0.04 


+U.27 
-0.13 


+0.28 
-0.15 



Table 2: Coefficients arising in the parametrization Hp with their different uncertainties. 

for practical purposes and for comparison with other parametrizations of the vector form factor, 
it is useful to give the first coefficients of the Taylor expansion 



m. 



m 



" dt 



m 



6 d^U 



" dt^ 



t=o 



t=0 



t=0 



A+, 

+ 2mlH'iO) = Al + (5.79+^;^^) x 10"^ , 

A^ + 6mlH'{0)A^ + 3m^if"(0) = A=^ + 3 (5.79^^;^^) x 10"^+ 
+ (2.99+°:^?) X 10^^ . (A.4) 



The third order Taylor expansion is very accurate in this case, too. The maximal error with 
respect to the exact dispersive representation is 6%. 



A.3 Calculation of / 



K 



In order to extract |/+ (0)14^1 from the measurement of the K^^ decay rate 

r^+/o = Afj,Ho Sew (1+2A^*%,) \ff'\o)V^f 4+/o , Afx+fo = C^+/o m^+/o/(1927r=^) 



(A.5) 

one has to evaluate the phase space integrals /^+/o defined in terms of the scalar and vector 
form factors: 



(A.6) 
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Using the dispersive parametrization for the form factors, see Eqs. (13.1113.71) . it is possible 
for practical purpose, to approximate the phase space integrals, Eq. (IA.6I) . by a polynomial 
expansion in terms of the two parameters InC and A+ entering this parametrization. One obtains 

Ik+/o = Co + CiA+ + C2A+ + CsAl + CiAi + c^lnC + CqIuC^ + cjlnC^ + cglnC^ , (A.7) 



where the polynomial coefficients for the 4 phase- space integrals are collected in the following 
table: 





Co 


Cl C2 C3 C4 


C5 C6 C7 C8 




0.14126 
0.09061 


0.48960 1.35655 3.12372 6.14597 
0.29599 0.95764 2.36194 4.81855 


0.01724 0.00515 0.00120 0.00023 




0.14530 
0.09324 


0.53899 1.59964 3.94735 8.32532 
0.32606 1.13005 2.98682 6.53170 


0.01798 0.00545 0.00129 0.00025 



Table 3: Coefficients ofA^ and \nC, Eq. t\A.7\) in the polynomial expansion of the phase space 
integrals, Eq. ({A.6\l . 
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